Macroscopic Quantum Self-trapping and Atomic Tunneling in Two-species 

Bose-Einstein Condensates 



O 

o 
o 

(N 



o 



«4H 

o 



o 
o 



> 
o 

m 
o 
o 
o 



I 

O 

o 



X 
J3 



Le-Man Kuang^'^, Zhong-Wen Ouyang^ 
^CCAST (World Laboratory), P.O. Box 8730, Beijmg 100080, China 
^Department of Physics, Hunan Normal University, Changsha 4IOO8I, China* 

(February 1, 2008) 

We present a new theoretical treatment of macroscopic quantum self-trapping (MQST) and 
quantum coherent atomic tunneling in a zero-temperature two-species Bose-Einstein condensate 
system in the presence of the nonlinear self-interaction of each species, the interspecies nonlinear 
interaction, and the Josephson-like tunneling interaction. It is shown that the nonlinear interactions 
can dramatically affect the MQST and the atomic tunneling, and lead to the collapses and revivals 
(CR) of population imbalance between the two condensates. The competing effects between the 
self-interaction of each species and the interspecies interaction can lead to the quenching of the 
MQST and the suppression of the CR and the Shapiro-like steps of the atomic tunneling current. 
It is revealed that the interatomic nonlinear interactions can induce the coherent atomic tunneling 
between two condensates even though there does not exist the interspecies Josephson-like tunneling 
coupling. 

PACS number(s): Q3.75.Fi, 74.50.-Hr, 05.30. Jp, 32.80.Pj 



I. INTRODUCTION 

Recently, much attention has been paid to the inves- 
tigations of systems consisting of two weakly interacting 
Bose-Einstein condensates [1-17] due to the appearance 
of quantum interference [18-34] and new macroscopic 
quantum phenomena [35-37]. In principle, such conden- 
sate systems can be produced in a double trap with two 
condensates coupled by quantum tunneling and ground 
collisions, or in a system with two different magnetic sub- 
levels of an atom, in which case the two species conden- 
sates correspond two electronic states involved. The cou- 
phng between two condensates could be realized by the 
near-resonant dipole-dipole interaction. 

The first experiment [6] involving interactions between 
two condensates was performed with atoms evaporatively 
cooled in the \F = 2, Mf = 2) and |1, —1) spin states of 
^'^Rb. One of the latest experimental advances [3,7] in 
this direction is the realization of measurements of rela- 
tive phase in two-component Bose-Einstein condensates. 
In experiments at JIL A [3] , two condensates in two differ- 
ent internal atomic states are produced by using a single 
two-photon coupling pulse. The two condensates have a 
well defined relative phase. After a time during which 
the condensates evolve in the trapping potentials, the 
two condensates interfere through mixing coherently the 
two internal atomic states. Then, the relative phase of 
the two condensates is obtained from the spatial inter- 
ference pattern. The realization of measurements of the 
relative phase between two condensates opened the fasci- 
nating possibility of experimentally examing the phase- 
related phenomena in Bose condensates, such as atomic 
Josephson effect and macroscopic quantum self-trapping 
(MQST) [37]. 

Theoretical studies of such systems began in Ho and 
Shenoy's work [9] which shown that binary mixtures 



of condensates of alkali atoms have a great variety of 
ground state and vortex structures. Then, the stabil- 
ity and collective excitations of two-species condensate 
systems [1,2,4,5] have been extensively studied. More re- 
cently, Smerzi and coworkers [37] have shown that the in 
a system of two Bose condensates the quantum coherent 
atomic tunneling between two condensates induces two 
types of interesting effects. One is an atomic Josephson 
effect in Bose condensates, which [43] is a generaliza- 
tion of the sinusoidal Josephson effects familiar in su- 
perconductors. The other is macroscopic quantum self- 
trapping (MQST), which is a kind of a self-locked popu- 
lation imbalance between two Bose condensates. It arises 
because of the interatomic nonlinear self-interaction. The 
MQST has a quantum nature, involving the coherence of 
a macroscopic number of atoms in the two condensates. 
It has been known that the MQST depends upon the 
trap parameters, the total atoms and initial states of the 
system and is self-maintained in a closed conserved sys- 
tem without external drives. As pointed out in Ref.[37], 
it is easier to observe the MQST in Bose condensates 
than self-trapping phenomena in other systems, such as 
the single-electron Coulomb Blockade effect [38] arising 
from the Coulomb interaction between electrons, single 
polaron trapping in a medium [39] which arises from sin- 
gle electrons, interacting with a polarizable lattice, and 
external gravitational effects on He II baths [40,41]. 

In Ref.[37], Smerzi and coworkers only considered the 
MQST induced by interatomic nonlinear self-interaction 
in each condensate. However, in a system consisting of 
two Bose condensates, there are not only nonlinear self- 
interaction but also interspecies nonlinear interaction. 
Questions that naturally arise are, what is the effect of 
the interspecies nonlinear interaction on the MQST and 
quantum coherent atomic tunneling? Does interspecies 
nonlinear interaction strengthen or weaken the MQST 
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and the atomic tunneling current between them? 

In this paper, we present a theoretical treatment of 
the MQST and the quantum coherent atomic tunneling 
in a more general two-species Bose condensate system in 
terms of a two-mode approximate model and the rotat- 
ing wave approximation. Our treatment involves not only 
the interatomic nonlinear self-interaction in each species 
but also the interspecies nonlinear interaction. We find 
that the presence of the interspecies nonlinear interaction 
gives rise to new insight to the MQST and the atomic 
tunneling between the two condensates. This paper is 
organized as follows. In Sec. II, we establish our model 
and present an approximate analytic solution. In Sec. 
Ill, we discuss the collapse and revival (CR) phenomenon 
on population imbalance between two condensates. In 
Sec. IV, we investigate the MQST in the two-condensate 
system, and discuss the dependence of the MQST upon 
the initial states and the tunneling interaction and the 
nonlinear interactions. In Sec. V, we study quantum dy- 
namics of the atomic tunneling current and its dc char- 
acteristics. We shall conclude our paper with discussions 
and remarks in the last section. 

II. MODEL AND SOLUTION 

We consider a zero-temperature two-species Bose con- 
densate system in which the atoms interact via aa and 
bb and ab elastic collisions, and there is a Josephson-like 
coupling term denoted by a^b and ab^. In the formalism 
of the second quantization, Hamiltonian of such a system 
can be written as 

H = Hi + H2 + Hint + Hjos, (1) 

Hi=l rfxVit(x)[-|lv^ + F,(x) 

+[/iit(x)V;i(x)]Vii(x), (i = l,2), (2) 
Hint = C/12 I rfxV3I(x)^|(x)^i(x)4(x), (3) 

Hjos =aJ dx[t/>l(x)V)2(x) + V'i(x)V'|(x)], (4) 

where '(/'i(x) and ^/'|(x) are the atomic field operators 
which annihilate and create atoms at position x, respec- 
tively, they satisfy the commutation relation 

[Vi,(x),Vi](x')]=5i,<5(x-x'). (5) 

In Eq.(l), Hi and H2 describe the evolution of each con- 
densate in the absence of interspecies interaction. Hint 
describes interspecies collisions. Hjos is the Josephson- 
like tunneling coupling term. Atoms are confined in har- 
monic potentials Vi{x){i = 1,2) of frequencies Wj. In- 
teractions between atoms are described by a nonlinear 
self-interaction term Ui = Anh^af'^/m and a term that 



corresponds the nonlinear interaction between different 
condensates U12 = 47rS.^af2/m, where a|'^ is s-wave scat- 
tering length of condensate i and af 2 that between con- 
densate 1 and 2. For simplicity, throughout this paper 
we let ?i = 1 and assume that af^ = = a"^, and 
yi(x) = F2(x). 

It is well known that the above Hamiltonian can be 
reduced to two-mode boson Hamiltonian [30-32, 35,36] 
through expanding the atomic field operators over single- 
particle states [35]: 

^i(x) = Oi^iijv(x) -I- Vii(x), (6) 

where d| = / (ix0jAr(x)-(/'|(x) create particles with distri- 
butions <?f)ijv(x) with [ai,a|] = 1. The first term in the 
mode expansion (6) acts only on the condensate state 
vector, whereas the second term V'ilx) accounts for non- 
condensed atoms. Substituting the mode expansions of 
the atomic field operators into the Hamiltonian (1), re- 
taining only the first term representing the condensates, 
we arrive at the following two-mode approximate Hamil- 
tonian 

H = woialai + 03^2) + Qial^al + a^a^) 

+g{a\a2 + ajai) + 2xa\aia2a2, (7) 

where the frequency and the coupling constants are de- 
fined by 
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^o = pJ (ix[^|V0i;v(x)|^ + F(x)(|</.,^(x)|2], (8) 
q = Uol dx(|0i^(x)|2 + |02w(x)|2), (9) 
9 = ^ J c;x(0{^(x)<^2jv(x) + c^in{x)4^{x)), (10) 

X=lUi2 y dx|4^(x)|2|<^2iv(x)|^ (11) 

From Eqs.(6) and (7) we can see that the two- mode ap- 
proximation essentially consists in neglecting all modes 
except the condensate modes. At zero temperature, this 
amounts to ignoring the atoms which have left the con- 
densate mode due to the two-body interactions. In other 
words, what the two-mode approximation involves is only 
the first order effects of interactions. The mode expan- 
sion of the condensate function over single-particle states 
(6) makes the condensate shape not to be changed, this 
limits migration of condensed atoms from one condensate 
to the other. The constraint on the shapes of condensates 
implies that the two-mode approximation can be applied 
only for weak nonlinearity. The valid conditions of the 
two-mode approximation were demonstrated in Refs.[30, 
32,35], which indicate that this approximation provides 
a reasonably accurate picture for weak many-body in- 
teractions , i.e., for small number of condensed atoms. 
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For large condensates, the mode functions of conden- 
sates are altered due to the collisional interactions, and 
the two-mode approximation breaks down. As shown in 
Rcf. [30,32], a simple estimate shows that this happens 
when the number of atoms iVa*^ 3> tq, where a**^ is a 
typical scattering length and vq is a measure of the trap 
size. If we consider a large trap with the size tq = 10/im 
and the typical scattering length a®'^ = 5 nm, the two- 
mode approximation is applicable for N < 2000. This is 
the case which we consider here. We shall show that the 
MQST and atomic tunneling between the two conden- 
sates are strongly affected by the nonlinear many-body 
interactions. 

We note that the two-mode approximate Hamiltonian 
has the same form with that of a two-mode nonlinear 
optical directional coupler [42]. The two-mode Hamilto- 
nian (7) can not be exactly solved, but for weak nonlin- 
ear interactions a closed analj^tical solution can be ob- 
tained under the rotating wave approximation suggested 
by Alodjanc et al. [43]. 

In order to obtain an approximate analytic solution of 
the Hamiltonian (7) , we introduce a new pair of bosonic 
operators Ai and A2 by the following expressions:: 

V 2 \/2 

(12) 

where Ai and A2 are slowly varying operators, they sat- 
isfy the usual bosonic commutation relations: [Ai,^j] = 
0, and [Ai, A|] = 5ij with yi| being the hermitain conju- 
gation of Aj. Then the Hamiltonian (7) reduces to the 
following form 
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H = wN+ ^q[3N^ - {AIAi - AlA2f 



+g{A\A, - A\A2) + |iV2 

-xA\a^aIa2 + h', 



(13) 



where the detuning is given by w = ~ (x + the 
total number operator TV is a conserved constant which 
is given by 

N = a\ai + ala2 = A\Ai + ijia, (14) 
and H' is a nonresonant term which is given by 

H[ = lix- q){AfAle-'^^' + AfAle^^% (15) 

which oscillates at the frequency Ag. The account of the 
fast oscillating term results only in some addtional oscil- 
lations which play no essential role in the evolution of the 
measurable quantities specifying the macroscopic quan- 
tum phenomena of the two-condensate system, so that it 
is fully negligible. This means the rotating wave approx- 
imation [43]. After neglecting the nonresonant term H' , 
we get the following approximate Hamiltonian: 



HA=ujN + g{A\A^-A\A2) 

+ ^^q[ZN^-{A\A^-A\A2f] 

+ ]pcN^ - xA\AiAlA2. (16) 

In order to solve the Hamiltonian (16) wo introduce 
two Fock spaces of {Ai,A2) and {0,1,(12) in which the 
bases are defined by 



|n, m) = 



1 



il"irio,o). 



\/n!m! 
1 



\n,m) = ^^a\^a\^\Q,Q). 
V n!m! 



(17) 
(18) 



where n and m take non-negative integers. Obviously, 
Ha is diagonal in the Fock space of (^1,^2), and we 
find that 



HA\n,m) = E{n,m)\n,m), 
E{n, m) = u){n + m) + g{n — m) 



(19) 



+ ^(3g + x){'n + mf - ^q{n - mf 
—Xnm. (20) 



Consider two coherent states defined in Fock spaces of 
(^1,^12) and (01,02), respectively. 



|Q;i,a2) = -Doi(Q:i)£'a2(a2)|0, 0), 

|W1,W2)=-D4,(W1)£'1,(W2)|0,0), 



(21) 
(22) 



where Da,i{ai) and D^ {ui) are displacement operators 
defined by 



Dai{oii) = exp(ajdj + a*a\), 
(ui) = exp{uiAi + u*Al). 



(23) 
(24) 



Note the fact that |0,0) = |0, 0), we can find a useful 
relation to connect ]ai, 02) and ]mi, U2) with each other 



, , ,ai + a2 i(ai — a2) , 
I Q!i — ia2 ai + ia2 . 



(25) 
(26) 



Following the arguments of Bosc broken symmetry, we 
assume that the two condensates are initially in the co- 
herent states ]ai) and \a2), which are eigenstates of fii 
and 02, respectively. Then the wave function of the two 
species condensate system at time t can be explicitly ex- 
pressed as 



1$ 



00 ^ 

(i))=e-^^ ^ ^^^uUiu2) 



xe 



n,m— 
-iE(n,m)t 



-\/n!m! 
n, m). 



(27) 
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where 



ui = -^(ai + a2),U2 = -^(ai - a2) (28) 



N=\aA^ + \a2Y = \m? + \u2\\ 



(29) 



where we have used Eqs.(25), (26) in the derivation of 
Eq.(27). 



III. COLLAPSE AND REVIVALS OF 
POPULATION IMBALANCE 

In this section we show that the two condensate sys- 
tem under our consideration exhibits a collapse and re- 
vival phenomenon of population imbalance between two 
condensates. Denote the number difference of atoms be- 
tween the two condensates by 



D{t) = N^{t) - N2{t). 



(30) 



Then from Eq.(27) we can find that at time t, the num- 
ber of atoms in each condensate Ni(t) = (a|oi) is given 
by 



Ni{t) = ^{N - {-iy2\ui\\u2\cos[4:gt + e{t)] 

^g-2JVsin=l(,-x)t|^(^=1^2) 

where we have used the following symbols: 



(31) 



|ui le''^". ,0{t) = {(pu2 -Vui)+ U21 sm{q - x)t, (32) 



with U21 and being defined by 



U21 = \U2\ 



_i Jail sin(/9i T (-l)1a2| sin(^2 . 
= tan \- ; — — — ; }. 

|ai| cos(/3i T (-l)*|a2| COS952 

Then, the population difference is given by 
D{t) =2\ui\\u2\cos[4gt + e{t)] 



(33) 
(34) 



xe 



-2iVsin"(i(g-x)t 



(35) 



where N = jaip + |a2p is the total number of the 
atoms in the two condensates. Eq.(35) indicates that 
the population imbalance periodically oscillates with the 
time evolution. From Eq.(35) we can see that D{t) 
exhibits collapse and revival phenomenon which is a 
kind of nonclassical effect well known in the Jaynes- 
Cummings model [44] to describe interaction between a 
single-mode radiation field and a two-level atom. The 
CR is also found in a Bosc condensate system [35,36]. 
From Eq.(35) we see that the CR of the population im- 
balance in the two-condensate system depends on the 
tunneling interaction (g) and interatomic nonlinear in- 
teractions {q and x)- When g > |(? — x|/8, since the func- 
tion cos[4:gt + 6{t)] is the rapidly varying part in (35), so 



that the shape of the CR is determined by the envelope 
function exp[— 2iVsin^(i((7 — x)t]- The maximal revivals 
take place at time t = 2nw/\q — x\, where n is an integer. 
When g < |g-x|/8, the function cxp[-27Vsin^(i((7-x)i] 
becomes the rapidly varying part in (35), the CR then is 
determined by the envelope function cos[4:gt + 9{t)]. In 
Fig. 1 we plot the evolution of the population difference 
between the two condensates with respect to the time 
which is in units of |g' — x|, when the two condensates are 
in the initial state of ai = 5 and a2 = 4, and the tunnel- 
ing coupling is g = 25\q — x\- Fig- 1 clearly indicates the 
CR phenomenon of the population difference. 

It is worthwhile to note that when q = Xi the popula- 
tion difference (35) becomes 



D{t) = 2|wi!|u2|cos[4.gt+ ((/?„, - (p^J], 



(36) 



which is a simple sinusoidal oscillation, no CR occurs. 
The suppression of the CR can be explained by look- 
ing at the expression (35) which indicates that both the 
self-interaction (g) and the interspecies interaction (x) 
can induce the CR, but the CR produced by one can 
weaken that by another. It is the CR produced by the 
self-interaction completely counteracts the CR by the in- 
terspecies interaction that leads to the suppression of the 
CR. 

When there is no nonlinear interactions, i.e., q = 
X = 0, from Eq.(35) it is easy to find that the popu- 
lation imbalance has the same form with that of the case 
9 = X 7^ 0- This means that the CR vanishes when 
the nonlinearity vanishes. Hence, the CR of the popu- 
lation imbalance is a consequence of nonlinear interac- 
tions in condensates. The CR of the oscillatory trans- 
fer of atoms between the two condensates constitutes a 
novel macroscopic quantum phenomenon induced by in- 
teratomic nonlinear interactions for the two species con- 
densate system. 



IV. MACROSCOPIC QUANTUM 
SELF-TRAPPING 

In this section we are concerned with the MQST. The 

MQST effect is characterized by the nonzero time mean 
value of the fractional population imbalance between the 
two condensates defined by 



Pit) 



Ni{t)-N2{t) 
N 



(37) 



From Eqs.(35) and (37) we get that 

2|ui||U2| 



Pit) = 



N 



■ cos[Agt + 0{t)] 



xe 



(38) 



In order to investigate the MQST, we expand the above 
equation as 
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+ 00 



P[t) = '-e Jn{u2l)Im{N) 



N 



n,m— — oo 



X cos{[(n + m){q — x) + 45]f 

+ {fu2 - 'Pui)}, 



(39) 



where J„ (A) and /„ (A) are Bessel function and modified 
Bessel function. From Eq.(35) it is easy to find that when 
the tunneUng interaction and nonhnear interactions sat- 
isfy the condition: 



4:g = K{x-q), 



(40) 



where K is an integer, we can get a nonzero time- 
averaged value of population imbalance 



n=—oo 



N 



X C0S(<^U2 -fui), 



(41) 



which indicates the existence of the MQST. Eq.(40) is 
the condition under which the MQST happens. 

From Eq.(39) we can see that when the tunneling in- 
teraction vanishes and nonlinear self-interaction equals 
nonlinear interspecies interaction, i.e., g = and q = Xt 
we arrive at a constant population imbalance 



Pit) 



2|mi||«2| 
iV 



(42) 



This is a time-independent state, called the self-trapping 
stationary state, which is the consequence of competing 
bctwcicn nonlinear self-interaction and nonlinear inter- 
species interaction. 

When there exists the tunneling coupling, i.e., g 0, 
from Eq.(40) we can find the critical value of the tunnel- 
ing coupling at which the MQST happens gc = I? — x|/4. 
This critical value gc depends upon only the difference be- 
tween the nonlinear self-interaction and interspecies non- 
linear interaction, not the nonlinear self-interaction and 
interspecies nonlinear interaction themselves. Therefore, 
it becomes possible that the MQST occurs only when the 
tunneling coupling equals or exceeds the critical value gc- 

In Fig. 2 we plot the time evolution of the fractional 
population imbalance when the two condensates are in 
the initial state of ai — 10 and a2 = for (a) K = 1, 
(b) K — 20. Here the time is in units of — x|. From 
Fig. 2 (a) and (b) we can see that the weaker the tunnel- 
ing coupling {g), the more apparent the MQST becomes. 
This indicates that the MQST is an effect induced by 
interatomic nonlinear interactions q and x not the tun- 
neling interaction (g). In what follows we shall discuss 
the dependence of the MQST in detail upon the initial 
states and nonlinear interactions for specific cases. 



We here discuss the dependence of the MQST on the 
initial states of the two condensates in the following four 
cases. 

Case 1: \oii\ = \a2\, ^ '^2- In this case, the two 
condensates initially have the same number of atoms but 
diflferent phases. Prom Eq.(32) and (33) we get that 
= iVcos2[((^i - ^2)/2], |M2p = iVsin2[(^i - ^2)/2], 
and M21 = — iVcos(</?i — ^2)- Making use of Eq.(28), 
from Eqs.(35) and we can see that the fractional popu- 
lation periodically evolves with respect to time, p{t) ^ 
except that ^pi — ^2 = nn, where n is an integer. If 
^g/{x — q) = K (an integer), we get the locked popula- 
tion imbalance 



p = I sin((^i - (^2)1 cos(i^„i - </?«2)e" 



-N 



-t-00 



X Jn{-NcOs{^i-^2))lK-niN), (43) 

n= — oo 

which implies that the MQST does exist, even if the two 

condensates initially have the same number of the atoms 
provided that they have different initial phases. And less 
the total atomic number is, stronger the MQST becomes. 

Case 2: \ai\ ^ \a2\, fi = <fi2- In this case the two 
condensates initially have the same phases but the dif- 
ferent number of atoms. From Eq.(38) we see that when 
4.9/(x ^ q) = K (an integer), the MQST occurs with the 
following locked population imbalance; 



P = 



|Q;i|^-h|Q2p 
2N 



+ 00 



e-^ E Jn{-2\ai\a2\)lK-n{N). 



(44) 

Case 3: N = jaip, \a2\ = 0. In this case, the sys- 
tem starts with all atoms being in one condensate. Mak- 
ing use of Eq.(28), we get that ui = U2 = ai/y/2, and 
6{t) =0. The fractional population evolution is given by 



+00 



Pit) 



-N 



E IniN)cos{[n{q-x)+4g]t, (45) 



So that when the tunneling interaction and the nonlin- 
ear interactions satisfy the relation ^g/ix — l) = K (an 
integer), we can see the appearance of the MQST with 
the locked population imbalance 



p = e-''lKiN) 



(46) 



Case 4- \cti\ = \a2\ and ipi = ip2- In this case the two 
condensates initially have the same number of atoms and 
the same phases. Prom Eq.(28) we have U2 = 0. Making 
use of Eq. (38) we can sec that the oscillations of the pop- 
ulation imbalance vanish, i.e., p(t) = 0, and no MQST 
occurs. This is in agreement with the result in Ref.[37]. 



A. The initial state dependence 



B. The dependence on nonlinear interactions 
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Then, wc turn to the dependence of the MQST upon 
the tunneling couphng (g) and the nonhnear interactions 
between atoms, which are described by the parameters q 
and X corresponding to self-interactions and interspecies 
interactions, respectively. 

Case 1: g = 0, q ^ 0, and x 0- In this case there 
is no tunneling interaction, but there exists interatomic 
nonlinear interactions. From Eq.(39) we find that the 
fractional population imbalance becomes 

n,m=— oo 

X cos{[(n + m){q- x}]t + {^U2 - Vm)}, (47) 

which indicates that no MQST occurs if the coupling of 
interatomic self-interaction does not equals that of inter- 
species nonlinear interaction, i.e., q ^ X- However, when 
the self-interaction equals the interspecies interaction, 
i.e., 9 = X) we can observe the self-trapping stationary 
state with a constant population imbalance p{t) = p{0). 

Case 2: g 7^ 0, g x — 0. In this case, we consider 
only the effect of the tunneling coupling while the inter- 
atomic nonlinear interactions are not involved. Eq.(38) 
tells us that the MQST does not occur, although there 
exists oscillations of the population imbalance between 
the two condensates. This further confirms the validity 
of Smerzi et al.^s conclusion [37] which the MQST arises 
from the interatomic nonlinear interaction. 

Case 3: g ^ 0, and 9 = X 7^ 0- this case there exist 
both the tunneling interaction and the nonlinear interac- 
tions, but self-interaction equals interspecies interaction. 
Prme Eq.(38) we can find that the population imbalance 
exhibits a simple oscillation with 

Pit) = M!f^e-^ cos[Agt + (^„, - <p„J], (48) 

which means that the MQST vanishes. 

Case 4: g ^ 0, q ^ 0, x = 0, ot g ^ 0, q = 0, x ^ 0- In 
this case, there exist the tunneling interaction and one 
of the self-interaction and the interspecies interaction. It 
is easy to see that the fractional population imbalance 
Eq.(38) reduces to 

Ht) = ^^e-^ g J„(«2i)/™(A^) 

n,m— — 00 

X cos{[(n + m)K + ig^ + {ipu^ - (puj}, (49) 

where k — q or — x- Eq.(49) reflects the fact that when 
4:g/K = n + m = K (an integer), the MQST happens 
with the nonzero p given by Eq.(41). This implies that 
both the nonlinear self-interaction in each condensate 
and the interspecies nonlinear interactions contribute to 
the MQST. Since the values of K to determine p have 
opposite signs for the self-interaction (x) and the inter- 
species interaction (q), the MQST produced by the self- 
interaction can weaken that by the interspecies interac- 
tion. It is the competition between the MQST induced 



by the nonlinear self-interactions of each condensate and 
that by the interspecies nonlinear interaction that leads 
to the quenching of the MQST in the above case 3. 

V. THE COHERENT ATOMIC TUNNELING 
CURRENT 

In this section, we study quantum dynamics of the co- 
herent atomic tunneling current between two condensates 
and its da characteristics, and discuss the influence of the 
initial state of condensates and the tunneling interaction 
and the nonlinear interactions. The coherent atomic tun- 
neling current between the two condensates is deflned 
by I{t) = Ni(t) - N2it). Making use of Eq.(31), it is 
straightforward to get that 

7(i) = -2|«i||«2|{45sin(0(i)) 

+{q - X)[\ui\'' MiQ - x)t - 0{t)) 

+ |u2psin((g-x)i + ^(t))]}. (50) 

This indicates that the atomic tunneling current period- 
ically changes, atoms periodically transfer between the 
two condensates with the time evolution. 

In order to see do characteristic of the atomic tunnel- 
ing current, we expand the atomic tunneling current (50) 
as the following expression 

+00 

I{t) = -2\ui\\u2\e-^ E {^9Jn{u2i) 

n,m=— 00 

-(9- X)[|W'lP-^«-l(W2l) - \U2\'^Jn+l{U2l)]} 

xI^{N) sin{[(m + n){q-x) + 4sf]t 

+ (¥'u2 - (51) 

which implies that when the tunneling coupling, and non- 
linear couplings satisfy the condition: 

Ag = K{x-q), (52) 

we get the dc component of the atomic tunneling current 
with the following form, 

Idc{K) = -2|ui||u2|(g-x)sin(v?u2 - V>ui)e~'^ 

X E {'2KJn{u2l) -[\ui\^Jn-l{u2l) 
n— — 00 

-\U2\^Jn+l{U2l)]}lK-n{N), (53) 

where K is an integer. This indicates that the dc compo- 
nent of the atomic tunneling current exhibits a step struc- 
ture with respect to the integer K. This step structure 
is a resonant phenomenon among the tunneling interac- 
tion and nonlinear interactions with the resonant condi- 
tion given by Eq.(51). It is the analogue of the Shapiro 
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steps observed in the superconductor Josephson junction 
[45], so that we call the steps in the step structure of 
the dc component of the atomic tunneling current the 
Shapiro- like steps. In what follows wc discuss in detail 
the dependence of the atomic tunneling current and the 
Shapiro-like steps upon the initial states and nonlinear 
interactions for some specific cases. 

A. The initial state dependence 

In this subection we discuss the initial-state depen- 
dence of the atomic tunneling current and the Shapiro- 
like steps for the following four cases. 

Case 1: \a\\ = \a2\, ^pi ^ ^P2- In this case, the two 
condensates initially have the same number of atoms but 
different phases. From Eq.(50) we can get the expression 
of the atomic tunneling current 

I{t) = -N\ sin(^i - (^2)|{4gsin(^^(^)) 

+(9 - x)^[cos2 - (^2) sin(((7 - x)t - 0{t)) 

+ siv? ^((^1 - ^2) sin((g - x)t + e{t))\), (54) 

where 

0{t) = {ipm - - Ncos{ipi - </?2)sin((g - x)t- (55) 

Form Eq.(54) we can obtain the dc component of the 
tunneling current with the following result, 

Idc{K) = -N{q - x)\ sinifi - (f2)\sm(Lpu2 - LpuJ 
xe-^ ^ {2i^J„(-iVcos(^i-^2)) 

n= — 00 

-A''[C0S^ ^(V?! - (f2)Jn-li-N COs{(pi - (^2)) 

-sin^ i((^i - ip2)Jn+i{-Ncos{ipi - <f2))]} 

XlK-n{N), (56) 

where n is an integer. Eq.(56) gives rise to the Shapiro- 
like steps of the atomic tunneling current. 

It is interesting to note that when the initial phases 
of two condensates satisfy the condition: fi — ip2 = nn, 
where n is an integer, we can find zero atomic tunnel- 
ing current. This means that the blockade of the atomic 
tunneling happens. 

When the initial phases satisfy the relation: (pi — (p2 = 
(2n -|- 1)77/2, n being an integer, we find that 

I{t) = -'igNsm{(pu2 - <fui) 

+{q - x)A''cos((/3„2 - ¥'«i)sin(2(g - x)t), (57) 

which indicates that the atomic tunneling current is a 
simple superpositon of a alternating current with the si- 
nusoidal oscillations and a dc current. The dc component 
is 



Idc = -4:gNsm{(pu^ -Vui), (58) 

which implies that the dc component of the atomic tun- 
neling current depends only upon the tunneling coupling, 
it is independent of the nonlinear interactions in two 
condensates, and increases linearly with the tunneling 
strength (g) and the the total number of the atoms (A''). 
In this case no Shapiro-like steps appears. 

Case 2: \oii\ ^ \a2\, = <f2- In this case, the two 
condensates initially have the same phases but the dif- 
ferent number of atoms. The atomic tunneling current is 
given by 

I{t) = -\a\-al\{^gsin{e{t)) 

+ - X)[{ai + a2f sin((g - x)t ~ 0{t)) 

+{ai-a2fsm{q-x)t + em}, (59) 

where 

e{t) = -2aia2 sm{q - x)t. (60) 

And the dc component of the atomic tunneling current 
has the following form, 

Idc{K) = -\al- al\{q-x) sin(<^„2 - (^„J 
xe"^ {2KJn{-2aia2) 

n= — oo 

-i[(ai - a2)^J„_i(-2aia2) 

-{ai - a2fJn+i{-2aia2)]}lK-n{N). (61) 

In FIG. 3, we plot the time evolution of the atomic 
tunneling current between the two condensates. Results 
are shown for the case oi g = 0.25 and g — x = 0.1 when 
the two condensates are in the initial state with ai = 5 
and a2 = 4. FIG. 3 indicates that the atomic tunneling 
current exhibits complicated oscillating behaviors. 

Case 3: N = |ai|^, \a2\ = 0. In this case, the system 
starts with all atoms being in one condensate. Taking 
into account Eqs.(32) and (33), from Eq.(50) we find the 
atomic tunneling current to be 

Iit) = -iq^x)N^sm{<l-x)t, (62) 

which is a pure sinusoidal alternating atomic current with 
the period T = 2Tr/\q — x\ without dc component. It is 
worthwhile to note that the atomic tunneling current (62) 
depends on only the difference between the self-coupling 
q and the interspecies coupling x, it is independent of the 
tunneling coupling g at all. Thus, we can conclude that 
the nonlinear interactions can induce the atomic tunnel- 
ing, even if there is no tunneling coupling between two 
condensates. 

Case 4: \a\\ = \a2\ and = ip2- In this case, the 
two condensates initially have the same number of atoms 
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and phases, wc sec that the atomic tunneUng current 
vanishes, and no Shapiro-Uke step occurs. 

The above analyses indicate that the atomic tunnel- 
ing current and the Shapiro-like steps strongly depend 
on the initial number of atoms in each condensate and 
the initial phase difference between the two condensates. 

B. The dependence on nonlinear interactions 

In what follows we show that the interatomic nonlinear 
interactions also significantly affect the atomic tunneling 
current and the Shapiro-like steps. 

If we consider only the influence of the tunneling cou- 
pling while the interatomic nonlinear interactions are not 
involved, i.e., 5 ^ 0, g = x = 0, we find that 

I{t) = -8g\ui\\u2\sm[4:gt + {(fu2 - fm)]- (63) 

This atomic tunneling current is a pure sinusoidally alter- 
nating atomic current with the period T = 7r/2(/. without 
dc component. Hence no Shapiro-like steps appears. 

It is interesting to note that when the interatomic non- 
linear interactions are involved, but g = x ^ 0, we get 
the same results with those of the case of g = x = for 
the atomic tunneling current and the Shapiro-like steps. 
This indicates that the contributions of the nonlinear self- 
interaction {q) in each condensate to the atomic tunnel- 
ing can counteract that of the interspecies nonlinear in- 
teractions (x)- It is the competition between the nonlin- 
ear self-interaction and the interspecies nonlinear inter- 
action that leads to a simple form of the atomic tunneling 
current and the suppression of the Shapiro-like steps. 

In particular, we note that when g = 0, i.e., there is 
no tunneling coupling, and q find tli^t 

I{t) = ~2\ui\\u2\iq - x)[|ui|' sin((g - x)t - e{t)) 

+ \u2\^sm{{q-x)t + em, (64) 

where 9{t) is given by Eq.(32). In FIG. 4, we display the 
time evolution of the atomic tunneling current (64) when 
the tunneling coupling vanishes. Results are shown for 
the case of q — x = 0.1 when the two condensates are in 
the initial state with ax = 10 and a-^ = 5. From FIG. 4 
we see that the evolution of the atomic tunneling current 
exhibits the CR phenomenon. 

Eq.(64) implies that the atomic tunneling current is 
nonzero, but no the Shapiro-like steps appears. From 
Eq.(64) we can see that even though there is no the tun- 
neling coupling between the two condensates {g = Q), the 
atomic tunneling between the two condensates may hap- 
pen. This atomic tunneling is completely induced by the 
nonlinearity of interatomic interactions which are char- 
acterized by interatomic collisions [q and x)- Therefore, 
we may conclude that the nonlinearity of interatomic in- 
teractions in the two condensates can lead to the atomic 
tunneling between the two condensates. 



VI. CONCLUDING REMARKS 



We have studied the MQST and the quantum coherent 
atomic tunneling in a two species Bose condensate sys- 
tem in the presence of nonlinear self-interaction of each 
species, the interspecies nonlinear interaction, and the 
Josephson-like tunneling interaction, and have given new 
insight to the MQST and the atomic tunneling. We have 
shown that the interatomic nonlinear interactions in the 
two condensates induce not only the MQST but also the 
CR of the population difference between two condensates, 
The CR phenomenon can be considered as a novel macro- 
scopic quantum effect. We have indicated that the non- 
linear interactions significantly affect the atomic tunnel- 
ing, and the Shapiro-like steps of the atomic timneling 
current. Comparing with Smerzi and coworkers' work 
[37], the present work involves the interspecies nonlinear 
interaction. The involvement of the interspecies nonlin- 
ear interaction gives rise to new characteristics on the 
MQST and the atomic tunneling. We have shown that 
when both the nonlinear self-interaction (g) and the in- 
terspecies nonlinear interaction (x) present at the same 
time, the atomic tunneling dynamics and the MQST and 
the Shapiro-like steps depends upon the difference (g— x)> 
not q and x themselves. We have also found that the 
interspecies nonlinear interaction generates the MQST 
at the same level with nonlinear self-interaction. How- 
ever, contribution from the interspecies interaction to the 
MQST and that from the self-interaction weaken them- 
selves with each other. It is the competing effects be- 
tween the nonlinear self-interaction in each species and 
the interspecies nonlinear interaction that leads to the 
quenching of the MQST and the suppression of the CR 
and the Shapiro-like steps of the atomic tunneling cur- 
rent. Especially, we have revealed that the nonlinearity 
of interatomic interactions in the two condensates can 
induce the coherent atomic tunneling between two con- 
densates occurs, even though there does not exist the 
Josephson-like tunneling coupling. It should be men- 
tioned that these results are obtained under the two- 
mode approximation and the rotating wave approxima- 
tion, so they are valid for weak nonlinear interactions be- 
tween atoms. Finally, It should be noted that in order to 
observe macroscopic quantum phenomena such as MQST 
and Shapiro-like steps in Bose condensates, one has to 
control various interactions between atoms. This control- 
ling can be carried out through manipulating interatomic 
scattering lengthes. Several theoretical and experimen- 
tal approaches [46-51] to alter the scattering length have 
been proposed. In particular, recent experiments on Fes- 
hbach resonances in a Bose condensate [46,47] have indi- 
cated that the scattering length of ultracold atoms can be 
altered through Feshbach resonance. These experimental 
progresses provide the possibility to observe the MQST 
and Shapiro-like steps in Bose condensates. 
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Figure Captions 

FIG. 1. Diagram of the time evolution of the popula- 
tion difference between the two condensates. The time 
is in imits of I? — x|. Results are shown for the case of 
g = 25|g — x|, when the two condensates are in the initial 
state with ai = 5 and a2 = 4. 



FIG. 2. Diagram of the time evolution of the fractional 
population imbalance. The time is in units of \q — x\- 
Results are shown for (a.) K = 1, and (b) K = 20 when 
the two condensates are in the initial state with ai = 10 
and a2 = 0. 

FIG. 3. The atomic tunneling current between the two 

condensates as a function of time t (in arbitrary units). 
Results are shown for the case oi g = 0.25 and q — x = 0-1 
when the two condensates are in the initial state with 
ai = 5 and a.2 = 4. 

FIG. 4. The atomic tunneling current between the two 
condensates as a function of time (in arbitrary units) 
when there does not exist the tunneling coupling. Re- 
sults are shown for the case oi q — x = 0.1 when the 
two condensates are in the initial state with ai = 10 and 
Q!2 = 5. 
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